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Reconstructing events with missing transverse 
momentum at the LHC and its application to spin 

measurement 



In this article we discuss the measurement of spin at the LHC, in events with two unknown 
four-momenta. Central to this problem is the identification of spin-dependent kinematic vari- 
ables and the construction of a statistical test that can distinguish between different spin 
hypotheses. We propose a method for reconstructing kinematic variables that depend upon 
the unknown momenta. The method is based upon a probabilistic reconstruction of each event, 
given the masses of the final and intermediate states and the cross-section of the assumed hy- 
pothesis. We demonstrate that this method can distinguish between two spin hypotheses for 
a specific process, even after mass uncertainties and Standard Model backgrounds are taken 
into account. We compare our method with another that only utilises the observable momenta 
of each event. We will show that our method permits an improved discrimination between 
hypotheses, with a reduced probability of error. 

1 Introduction 

Since the end of March 2010 high-energy collisions have been taking place at the Large Hadron 
Collider (LHC), ushering in a new era of particle physics. Over the coming decades the 
experiments at the LHC will test the Standard Model (SM) in a new energy regime and will 
search for the first conclusive evidence for physics beyond it. One possible signal for this new 
physics, which is motivated by the existence of dark matter, are events with large missing 
transverse momentum. Such events may have the following topologjQ: 
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where constituents, //, of the cohiding protons interact to produce a pair of new particle 
states AA. These subsequently decay into an electrically neutral and colourless state x (the 
DM candidate), that will escape detection, and a set of visible SM states VV. Because the 
four-momenta of the two x particles are unknown, the kinematics of these events cannot be 
reconstructed. This loss of information will make it difficult to measure the masses, spins 
and couplings of these particle states. Such measurements, however, will be crucial if the 
physics underlying these signals is to be understood. In this paper we focus on the issue of 
spin measurement in events of this topology. 

In order to determine the spin of either A or one must measure the probability density 
of some kinematic quantity. This distribution will then be compared to those predicted by 
several hypotheses, e.g. for different spin assignments to the particles in the event, and one 
must construct some statistical test in order to decide which hypothesis is correct. Central 
to this problem, therefore, is the identification of kinematic quantities, and an associated 
statistical test, that will permit a spin determination with a small probability of error. The 
ultimate goal is to identify an 'optimal' set, which will give the smallest probability of error 
for a fixed luminosity. 

Previous studies of spin determination at the LHC have only considered kinematic quanti- 
ties that are observable — i.e. those that can be uniquely determined for each event. Examples 
include Lorentz invariants formed from the visible momenta of decay products [IHS], and 
differences in pseudo-rapidities [5]. These variables, however, may not be optimal. 

Alternative variables that depend upon the spin are the polar and azimuthal angles of 
production and decay, measured relative to some physical axis in the centre of mass frame and 
the rest frames of decaying particles. The energy dependence of the cross-section is another 
spin-dependent variable. It has been shown that using these variables one can determine the 
spin uniquely, in certain processes [IH^ . The difficulty with these variables is that they depend 
upon the unknown final state momenta, and so cannot be reconstructed event- by- event. This 
does not imply, however, that one cannot obtain some information regarding the probability 
density of these variables. 

Consider some variable, O, with probability density p{0\H) for some given hypothesis, 
H. (In principle H must specify all relevant physical quantities, including spins, masses and 
couplings.) It follows that: 

p{o\h) = J2p{o\v,h)p{v\h) (2) 

V 

where the sum is over each point, in the phase space of the visible final state. Therefore, 
suppose we construct, for each hypothesis, the following distribution for a given set of N 
events: 

N 

qH{0) = Y,p{0\E,,H) (3) 

1=1 

where the sum is over each of the events Ei in the set. The conditional probabilities p{0\Ei, H) 
include all of the information that we have about an event, including the observed visible 
momenta and any information regarding the masses of A and x- If H is the true hypothesis, 
this distribution will have the following form in the large statistics limit: 

lim qniO) ^ y^p{0\V,H){Nv) (4) 

V 

where the sum over events has been rearranged into a sum over phase space points, with 
{Ny) the average number of events at the point V . Since {Ny) = Np{V\H), it follows from 
Eqns. ^ and (g]) that 

lim qH{0) ^ Np{0\H) (5) 
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Figure 1 : Defines the convention for momentum assignment in the process of Eq. ([Ij . 



In general, qn will tend to a different distribution if iJ is a false hypothesis. Hence, in principal, 
one can identify the true hypothesis by comparing the distribution qniO) with Np(0\H). It 
is this idea that we will explore in this paper. 

The remainder of this paper is organised as follows. In Sec. [5] we will demonstrate how one 
can construct, for a given hypothesis and kinematic variable, the distribution qniO). We will 
then discuss how to construct a statistical test to compare different hypotheses and decide 
which is true. 

In Sec. [3] we will apply our method to the discrimination of SUSY and UED models 
in slepton pair-production. For reference, we will compare our method with another that 
currently exists in the literature. We note that the purpose of this study is not to perform a 
fully rigorous analysis of how accurately our method will perform at the LHC. Such a study 
will require a full understanding of the detector and backgrounds, combined with a detailed 
study of systematic uncertainties. This is beyond the scope of this work. Rather, our goal is 
to perform a preliminary study, in a somewhat idealised case, in order to potentially identify 
the best method to use for a spin measurement. 

Finally, in Sec. |4] we will summarise our results and make some suggestions for further 
developments. 



2 The method 

This section is composed of two parts. In the first we will discuss how to construct the 
distribution qn for a set of events and a given hypothesis H. Then, in the second part, we will 
discuss how qn can be used in a statistical test to discriminate between different hypotheses. 

2.1 Constructing qh 

Let us consider a set of events with the topology of Eq. ([T]), where the initial and final state 
momenta are labelled as shown in Fig. [1] From this set we wish to construct the distribution 
<1h{0), defined in Eq. ([3]), for some kinematic quantity O and hypothesis H. To do this we 
perform the following procedure. For each event, Ei: 

• We determine the space of all possible four-momenta k and k that can be assigned to 
Ei, consistent with the hypothesis H. We will call this the space of solutions. 

• We construct the probability density p{k, k\Ei,H) on this space. 

• We integrate p{k,k\Ei, H) over the space of solutions to determine the distribution 
p{0\E,,H). 

We then sum the p{0\Ei, H) from each event to obtain qn- We will now discuss each of these 
points in more detail. 
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2.1.1 The space of solutions 

The final state of an event Ei is parameterised by tlie eight unknown components of the two 
four- vectors k and k. Not ah of these components are independent. Due to the conservation 
of four-momentum there exists the constraint: 

PA+PB = {E,0,0,E^) ^p + k+p + k (6) 

where E and Pz arc the (unknown) total energy and momentum, respectively, in the LAB 
frame. Hence the transverse components of k and k must satisfy: 

kr + kT = firp = -{pT +pt) (7) 

We will further assume that the masses and m^, of the states A and x respectively, 
have been measured. Several methods have been proposed to measure these masses in a model 
independent manner, and we refer the reader to the relevant literature jlO| . Hence, knowing 
the mass, one can apply the on-shell constraint to the unknown momenta: 

fc2 = ml F ^ ml (8) 

The intermediate state A will have some finite width, r,^. It is unlikely that a reliable 
measure of Tj[ will be possible at the LHC. Given this we will make the narrow width approx- 
imation for this state, which has the additional advantage of providing two further on-shell 
constraints: 

(p + fc)2=m3^ {p + k)^^ml^ (9) 

Provided T^/uia <C 1, this is likely to be a good approximation. 

The six constraints given by Eqns. ([7]), ([5]) and (jH]) reduce the number of free parameters 
in the final state to two. Without loss of generality we choose these free parameters to be the 
transverse components fc^. By solving the constraint equations one finds, for a given kx-, two 
solutions for the four- momenta k: 

fc(z) = (fc°W,fcT,fcL(l)) (10) 



k%i) = Jml + kT-kT + kl{i) (11) 



where i G {0,1}. We have defined 5w? = — m'^{p,k), where rriTia^b) is the transverse 
mass of the four- vectors a and b: 

m|(a, 6) + 6^ + 2(a?,6?, - AT • 6t) (13) 



and aj^ ~ \/ -\- olt ■ a.T is the transverse energy. Through the substitutions p ^ p and k ^ k, 

with kx = frp — kx^ one obtains the equivalent solutions for k. For each choice of kx there are 
four discrete solutions, corresponding to the two different sign choices for k^ and k^. 

A real solution exists only if the two following constraints are satisfied: ■m'^{p,k) < 
and rn^{p,k) < m^. By choosing an orthogonal set of co-ordinate axes in the transverse 
plane such that px = |pt|(1, 0), the first constraint can be written in the form: 

p2fc2 + (pO )2^2 _ 2ci\pT\k:, - Cq < (14) 

where the constants cq and ci are given by: 

co^cl- ip'^xfrnl (16) 
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whilst the second has the form: 

{Prfkr ■ kr - {pr ■ krf + 2di • fcr + rfo < (17) 

where 

di = (ci + PT ■ ^t)pt - (Pt)¥t (18) 
do = {Prf^T ■ ^ (PT ■ frf ^ 2cipT ■ - cq (19) 

Thus kr is constrained to he within the overlap region of two ellipses, which are defined by 
Eqns. ([T4| and ([T7|) . Thus the space of solutions associated with an event Ei is the set of four 
discrete solutions associated with each kr that lies within this overlap region. 



2.1.2 Constructing a probability measure on the space of solutions 

The probability measure p(fc, k\Ei, H) is given by the differential cross-section: 



da = dx dx'^fa{x)h{x')dcrab 
JO Jo „ I. 



(20) 



where fa are the proton's parton density functions (PDFs) and the sum is performed over all 
parton flavours a. We define x and x' such that: 

PA = Ebx{l,Q,Q,l) (21) 
PB=Ebx'{l,Q,Q,-l) (22) 

where E^ is the beam energy. For the LHC we will assume the full design beam energy of 7 TeV. 
The cross-section d& for the hard scatter is given by: 

^^-^ = 2^ ( n T^^O^ ) i'^^fSHpA+PB J2 (23) 

where -^/s is the total energy in the CM frame and a general four-momentum in the final state 
is written a;. The matrix element for the process, which is specified for a given hypothesis H, 

is Mab- 

Making the narrow width approximation for the intermediate state A forces A4ab to have 
the following structure: 

IXabP = \Mab\''S{{p+k)^ -m^^)S{{p+kf -ml^) (24) 

where ah is a dimensionless quantity. Inserting Eq. (j24[) into Eq. (|23p and integrating over 
the components of k and kL, one finds, up to irrelevant constant factors: 

dV ^30 \Mab\^ d^kr 1 (2). ^ sr \ 

'^""^ °^ mm ^ ■ /fc^fc^|sinh(Ay)sinh(Ay)|^ ^"'^ ^'^^ 

where the sum over k^ and k^ corresponds to a sum over the solutions given in Eq. (1121) . The 
rapidity differences Ay and Ay are given by: 

1 (^+PL k^_kL\ 1, ff+PL P-fcL\ 
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Hence the total differential cross-section has the form (again up to constant factors) : 



^ s ■ fcOfcO |sinh(Ay)sinh(Ay)| ^ ^ 

where x and x' are fixed by the requirement: 

Eb{x + x\0,0,x - x') ^p + k+p + k (28) 
0<.t,.t'<1 (29) 

Hence the probability measure p{kT, i, j\E, H) is given by: 

P^'^-^'^^^^^''^'' ^--^ ^> fcg.4|sinh(A,)sinh(A^)| ^''^^ 

where ?, j G {0, 1} correspond to the different k^ and k^ solutions at each point. The constant 
N is defined such that Ylii j j p{kT,h3\E,H)(PkT = 1- 

2.1.3 Constructing p(0|£;,i7) 

Given the space of solutions defined in Sec. 12.1. H and the probability measure of Eq. ([50]) . one 
can obtain the probability density for any kinematic quantity O through a straightforward 
change of variables. In practice we simply integrate p(k,k\E,H) over the region of solution 
space for which O takes a given value: 

p{0'\E,H)^ J dkTj2p{kT,t,j\E,H)diOikT)-0') (31) 
where the integral and sum is performed over the entire space of solutions. 



2.2 The statistical test 

The previous sections have detailed how to construct qn for a set of N events. Given this 
distribution for each hypothesis, one would then like to decide which hypothesis is true. For 
this, we must construct a statistical test. 

Since qn depends upon the random events in the set, qn is itself random and so will 
be distributed according to some probability density. For events distributed according to the 
hypothesis H let us denote the probability density of qH as f{qH\H). (That is, f{qH\H) is the 
probability density for qn in events generated and reconstructed with the same hypothesis. In 
statistics terminology, one can also consider f{qH\H) to be a likelihood.) Let us assume that 
there arc / hypotheses Hi, with i £ {1, • • • J}, each with an associated probability density 
f{qHi \Hi). If we knew each of these probability densities we could construct a statistical test 
as follows. 

Suppose that we are given a set of N events, and we are asked to determine which of the 
/ hypotheses is true. For each hypothesis we could calculate gjj. , as detailed above, and then 
calculate the associated probability, given the hypothesis, fi = f{qHi\Hi). We then reject the 
hypothesis Hj in favour of Hi if 

rij = — > Cij (32) 
Jj 

where the constants Cij are chosen so as to minimise errors (i.e. rejecting Hi when Hi is true, 
or accepting Hi when Hi is false). We will call r^j the likelihood ratio. 

Unfortunately we cannot calculate the probability densities f{qH\H) analytically. Instead, 
they will be determined numerically, through Monte Carlo simulation. To simplify matters 
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Sa 


R{A) 


X 






SUSY 
UED 





1/2 


(0,2,1/2) 
(0,2,1/2) 




1/2 
1 


(0,0,0) 
(0,0,0) 



Table 1: Table of spin S and representation R assignments in the SUSY and UED models. (Note 
that in UED ^ is a Dirac field, with both left and right handed components transforming in the same 
representation.) 



we will divide O into AI discrete bins, where M is made sufficiently small so as to simplify the 
calculation, without losing significant detail in the shape of qh- Thus for each set of events 
one calculates the M discrete bin values q}^, with i G {!,••• ,M}. We will denote this set of 
M values by the vector Qh = {q^, ■ ■ ■ , Qh}- Hence the problem of determining / is reduced 
to finding the Af -dimensional probability distribution f{qH\H). 
Since q\j is the sum over many random variables: 

N 

g}, = ^p(0'|i?„i/) (33) 

one expects from the central limit theorem that it will have a Gaussian distribution. (This 
will hold provided iV ^ 1 . For any process that may be viably studied at the LHC one might 
expect TV ~ 0(10^), so this assumption seems reasonable.) It seems sensible therefore to 
assume an i\/-dimcnsional Gaussian distribution for /(qnlH): 



fiqH\H) = ^ exp ( - liqn - flHfV^'iqH - m)) (34) 

where is the average value in each bin and Vh is the M x M correlation matrix. In general, 
since qn is a smooth distribution, there will be significant correlations between the bin values. 
The ^M{3 + M) parameters of this Gaussian distribution will be estimated using Monte Carlo 
simulation. We note, however, that it follows from the definition of qn that: 

/i' = (9^) = {N)piO^\H) (35) 

which provides an important cross-check of our method. 

2.3 Summary 

In this section we have discussed how one can construct qn^O) and then use this in a statistical 
test to discriminate between hypotheses. In order to illustrate the method we apply it to an 
example in the following section. 



3 A test of the method: Discriminating between models 
with different spin 

As a first test of the method we have attempted to discriminate between two models in the 
following process: 

qq^AA^ +X + + X (36) 

The models differ in their spin assignments and are based upon: supersymmctry (SUSY), 
where A is the scalar superpartncr of the muon and % is a Majorana fermion, and universal 
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(a) SUSY 



(b) UED 



Figure 2: A plot of the cos6'^ — cos0^ distribution in the SUSY and UED models. 

extra dimensions (UED), where A and x ^-re the first Kaluza-Klein exeitation of the muon 
and B gauge boson, respectively. The spin 5* and representation R under the SM gauge group 
SU{3) X SU{2) X U{1)y of A and x are summarised for eaeh model in Table [1] 

Several previous works have shown that a (partial) measurement of the spin of A is feasible 
in this proeess, both at a linear collider and at the LHC [6l|8l|9]. Let us summarise, briefly, 
how this can be achieved. 

The qq initial state of this process is forced to have a total angular momentum = ±1. 
(This is a consequence of the antiparallel momenta of the q and q together with their coupling 
to an s-channel vector boson.) Hence, in the SUSY model, the scalar pair (AA) must be 
produced in a P-wave state, in order to conserve angular momentum. These states carry 
orbital angular momentum and have a non-trivial angular distribution. Thus the SUSY matrix 
element has the form: 



where 9j{ and /3_4 are the polar angle of the A momentum, measured relative to the beam- 
line, and the speed of A, as observed in the CM frame. The factor k is dimensionless and 
independent of 9j{. This matrix element vanishes in the limits sin^.^ — > and Pj[ — >■ as 
a straightforward consequence of angular momentum conservation. If, however, A carries 
intrinsic angular momentum through a non-zero spin, total angular momentum can be con- 
served even with vanishing orbital angular momentum. Hence the (^^) can be produced in 
an S'-wave state, which has a matrix element that is non-zero for all 0yi and f3^. The precise 
form of the matrix element depends upon the spin and couplings of A. If A is spin-1/2 with 
vector-like couplings, as in UED, this matrix element has the form: 



Thus the distributions of both cos 9^ and are sensitive to the spin and can be used to 
discriminate between the SUSY and UED models. A more general analysis has shown that for 
a process of this kind, these distributions are sufficient to discriminate between the hypotheses 
Sa^O and 5^ > ig. 

Other variables are also sensitive to the spin, such as the azimuthal and polar angles of 
the fj,~ (/x"*") momenta in the rest frame of the A (A). Of particular interest are the following: 

• The azimuthal angle (j) between the decay planes of the A and A, defined in the CM 
frame. As pointed out in [S], has a distribution of the form: 



1-^ I SUSY = ksusy/33v sin^ 



(37) 



l-A^luED = '«UED(2-/33isin2 9a) 



(38) 




(39) 
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Hence one can infer tlie spin oiA from tlie higliest cos(2i0) component of the distribution. 

• The polar angles 0^ (^^) of the fi' (m^), measured in the rest frame of the A {A). These 
are measured relative to the direction of the boost from the LAB to the CM frame, as 
viewed in the appropriate rest frame. In the rest frame of A only its polarisation can 
break rotational invariance. Hence the observation of a non-trivial distribution in 0^ 
and 9^ indicates that A must be polarised and hence have a non-zero spin. For the UED 
model there is a large correlatiorH between these angles, whilst the SUSY distribution 
is flat, as shown in Fig. [2] 

Unfortunately, all of these variables depend upon the unknown x momenta, so they cannot 
be uniquely determined for an event at the LHC. However, we propose to employ our method 
to reconstruct their distributions and so discriminate between models. 

An alternative method has been proposed for use at the LHC, which depends upon the 
variable cos Ou [5] : 

cos0„ =tanh (^^^^ (40) 

where Arj^ is the pseudorapidity difference between the /i+ and fi" . This variable depends 

only upon the observable momenta and so is straightforward to reconstruct. Because it de- 
pends upon a pseudorapidity difference, it is invariant under longitudinal boosts and so is 
sensitive to the direction of the /i^ system in the CM frame. Since these states tend to be 
boosted in the direction of their parent A and A this variable is correlated with the polar 
angle 0^, and so is sensitive to the spin. We propose to use this method as a benchmark for 
our own. 



3.1 Simulating events at the LHC - signal, background and event 
selection 

The LHC experiments will search for events from this process using the signature of two 
opposite sign muons, in association with missing transverse energy. There will, however, be 
significant background sources of such events arising from SM processes. 

A realistic test of our method must take into account the effect of these backgrounds. 
However, a complete background analysis, combined with detector effects and associated 
systematic errors is far beyond the scope of this work. Instead we will perform our test only 
at parton level, without hadronisation and detector simulation. Furthermore, we will only 
include a limited number of backgrounds in this analysis. Previous studies [6] have identified 
several important sources, the most important of which is the irreducible SM background: 

qq^ 11^ + /.t" + + (41) 

This is the only one that we will take into accounlH. Despite these simplifications this test 
should still establish whether, in principal, our method improves upon the use of the cos On 
distribution. 

In order to isolate the signal events from background we follow the work of Barr and 
apply the following selection criteria: 

■^This correlation arises because the decay vertex involves a chiral coupling and there is a net polarisation of the 
initial state along the direction of the boost, B, from the CM frame to the LAB frame. This polarisation arises 
as follows. Because .4 is a doublet under SU{2) the process qq AA involves a significant contribution from 
s-channel Z exchange. Since the Z has chiral couplings to the quarks, this results in an overall polarisation of the 
initial state along the direction defined by the incoming quark's momentum. However, owing to the structure of 
the pdfs and the mass of A, this direction is correlated with the direction of the boost B. 

•^Using the cuts listed below, it has been found [6] that this process accounts for over 80% of the background 
signal. 
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• The transverse momenta of the muons must satisfy pTi > 40 GcV and pT2 > 30 GeV, 
where, by definition, pxi > pt2- 

• The invariant mass of the muon pair must satisfy: m^p > 150 GeV 

• The missing transverse energy must satisfy > 100 GcV 

• Mt2{M = 0) > 100 GeV 

where the kinematic variable Mt2 [IIIIH] is defined as: 

Mt2{M) = min ^max (mT(p, k), rriTip, fc))^ (42) 

where the minimisation is subject to the constraints k^ ^ k'^ = AP and kT + kT = frp- These 
cuts have been developed to eliminate most of the events from the SM di-boson processes: 
W~^W~ and ZZ production, but have not been tuned to maximise the ratio of signal to 
background. In addition to these criteria we also demand that the final state muons have a 
pseudorapidity |?7| < 2.5, which is required to ensure a high reconstruction efficiency in both 
ATLAS and CMS mUl]. 

Using MadGraph/MadEvent [TS] wc find that the leading order cross-section for the SM 
process given in Eq. (j41[) is, after applying the above cuts, 

asM = 0.5 fb (43) 

The signal events for both the SUSY and UED model were also generated using Mad- 
Graph/MadEvent. The UED model was implemented as a user-defined model, with particle 
widths calculated using the program BRIDGE [TB]. For this test we fixed the masses at 
mj\ ~ 350 GcV and — 50 GeV. For the SUSY model we found a leading order cross- 
section, after cuts, of 

^TsusY - 1.0 fb. (44) 

whilst the cross-section for the UED model takes a significantly larger value. Because the 
couplings of the new particle states cannot be known prior to a spin determination, the cross- 
section cannot be used as true spin (or model) discriminant. (Our method only relies upon the 
shape of the qn distributions to perform the discrimination.) Thus we simulate both SUSY 
and UED models by assuming the same cross-section for each and, to be conservative, we 
choose the smaller cross-section of the SUSY model. 

3.2 Implementing and testing our method 

For this test we performed a total of 2000 pseudo-experiments, each with a total integrated 
luminosity C = 100 fb^^, for each model. In each pseudo-experiment the number of signal and 
background events were selected at random from Poisson distributions with means (Ns) = 100 
and (Nb) = 50, respectively. The corresponding number of signal and background events were 
then generated and reconstructed using our method. 

When calculating qn for each hypothesis we made several simplifications in the probability 
measure p{k,k\E^ H), given in Eq. pO|) . Rather than use the full matrix element for either 
hypothesis, we instead used the matrix element for the production process only, which are 
given in Eqns. (|37| and ((38)) . The k factors were treated as dimensionless constants and 
were absorbed into the normalisation factor A/". In addition, we only include the PDFs for 
the up- type quarkfl which were calculated using the NNLO central MSTW set [17], at a 
scale Q = m^. These approximations were motivated by the desire to lessen the dependence 
of the reconstruction on unknown, mode-dependent factors, such as the structure of decay 

''We have found that there is no significant difference between using the PDFs for either the up- type or down-type 
quarks. 
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vertices and the strengths of couphngs. With these approximations we are effectively left with 
with a probabihty measure p(k, k\E, H) that is independent of these unknown factors. The 
approximations neglect some very mild dependence on ^/s that is actually present in the k 
factors, and spin correlations in the decays of A are also neglected. 

Because the number of signal and background events expected to pass the selection cri- 
teria are comparable, the reconstructed qn distribution will contain a significant background 
contamination. A background subtraction was performed in order to isolate the signal con- 
tribution to qn- To do this we determined the average background contribution, (g//), which 
was then subtracted in each pseudo-experiment to obtain a new distribution: 

Qh ^Qh - {qh) (45) 

These subtracted distributions were then used to perform the statistical tests discussed in 
Sec. 12.21 In particular, they were used to determine the means and the covariance matrix for 
the probability densities /(gsusvlSUSY) and /(quedIUED), using unbiased estimators. These 
probability densities were then used to calculate, for each pseudo-experiment, the likelihood 

/(IsusyISUSY) 



/(guEDlUED) 

3.3 Results 



(46) 



To begin, we will discuss our results from an attempt to reconstruct the cos 9^ distribution, 
which we divided into AI = 10 bins. To perform the reconstruction we assumed masses 
m-^ = and uiji = riiA, where rh^ = 50 GeV and to^ = 350 GcV are the true (pole) 
masses of x and A, respectively, that were used to generate the events. In Sec. 13.3.1] we will 
discuss the effect of errors in these assumed masses. 

In Fig. [3(a)] and Fig. [3(b)] we show the average g//(cos6'^) distribution for H = SUSY and 



H = UED, respectively. In each plot we show the distribution obtained from events generated 
for the SUSY model (filled circles) and UED model (open circles). The histogram in each plot 
shows the {Ns)p{cos9j^\H) distribution. Also shown in each plot are error bars, where the 
error cr* in the i"^ bin is given by the corresponding estimated variance: 



cr' ^ JV^ (47) 

It can be seen from these figures that events reconstructed with the true hypothesis have a 
qu distribution that is in agreement with p(cos 9ji\H). This is as we expected, and provides an 
important validation of our method. However, we observe that there is a small discrepancy for 
events from the UED model. We expect that this is caused by our neglect of spin-correlations 
in the reconstruction, as discussed in Sec. 13.21 

Crucially, we observe that there is a significant discrepancy between events from the SUSY 
and UED models when reconstructed with the same hypothesis. This demonstrates that the 
method can discriminate between the two models, at least in principal. In Fig. |4] wc show 
a plot of the r distribution for events from SUSY (hatched histogram) and UED (empty 
histogram). The SUSY and UED distributions are well separated, thus indicating that r is a 
potentially useful test statistic. To define our statistical test we choose our acceptance region, 
in which we accept the SUSY hypothesis instead of UED, as: 

ln(r) > 2 (48) 

From the distribution of Fig. 0] this test is expected to have a significancifl: 

a= / p(ln(r)|SUSY)dln(r) = 3.3% (49) 



^The significance of a hypothesis test is defined as the probability of an error of the first kind, i.e. the probability 
of rejecting the hypothesis of SUSY, if SUSY is the true hypothesis [15] . 
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Figure 3: The reconstructed qnicos 9j\) distributions for events from the SUSY (filled circles) and UED 
models (open circles). 
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Figure 4: The r distribution for events from the SUSY (hatehed histogram) and UED (empty histogram) 
models, eonstructed using the qnicosO^) distribution. The dashed hue and arrow indicate the acceptance 
region for the SUSY hypothesis, given in Eq. (|48| . 



and powei|f|: 



1-/3 = 1 



p(hi(r)|UED)dln(r) = 98.3% 



(50) 



To put this resuh into context we compared this with a test based upon the cos 6u distri- 



bution. In Fig. 5(a) we show the average cos On distribution for our sample of 2000 pseudo- 
experiments, for both SUSY and UED models, after the average SM background contribution 
has been subtracted. 

We constructed a statistical test for this method as follows. We assumed, for each hy- 
pothesis H, that the number of events nf in the i'^ bin, after subtraction of the averaged 
background, is a Gaussian random variable with a mean fif and variance [cr^Y that were es- 
timated from the set of pseudo-experiments. Using these Gaussian distributions we calculated 
for each pseudo-experiment the likelihood ratio 



rii 



p(rti, ...,r7,io|SUSY) 
p(ni, ...,nio|UED) 



n 

i=l 



,UED 



.SUSY 



exp 



SUSY 



UED 



2(^SUSY)2 



2(^yED)2 



(51) 

(52) 



The distribution of m is shown for both SUSY and UED models in Fig. [5(b)| The SUSY and 
UED distributions are separated, however there is a significant overlap between the two. If 
we define the acceptance region as: 

Hni) > (53) 

we find a statistical test with acceptance a = 14.8% and power 1 — /3 = 84.9%. Comparing 
this with Eqns. and ([50]) we find that our method reduces the probability of error by a 
factor of ^ 5. 



''The power of a hypothesis test to discriminate between the hypotheses is given by (1 — /3), where /3 is the 
probability of an error of tlie second kind, i.e. the probability of accepting the hypothesis of SUSY when SUSY is 



false [H]. 
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Figure 5: The cosOu and ru distributions for events from tiie SUSY (filled circles, hatched histogram) 
and UED (open circles and empty histogram) models. The dashed line and arrow indicate the acceptance 
region for the SUSY hypothesis, given in Eq. ([SS]). 
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(a) H = SUSY 
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(a) H = SUSY 
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Figure 7: The average joint (7i/(cos6'^, distribution for UED events, after background subtraction. 
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log(r) 

Figure 8: The r distribution for SUSY (hatched) and UED (empty) events, constructed from the joint 
(7h(cos 0^, /3^) distribution. The dashed line and arrow indicate the acceptance region for the SUSY 
hypothesis, given in Eq. (|54|) . 

Following this initial test we attempted to improve our model discrimination by including 
additional kinematic variables. In Fig. [S] we show the joint distribution qh {cos 6 [A) , for 
the SUSY events, reconstructed with both hypotheses after the average background has been 
subtracted. The distribution for the UED events is shown in Fig. [7l One interesting feature 
of these distributions is that the dependence is insensitive to the hypothesis H used in 
the reconstruction. As can be observed in Fig. [6] the distribution for the SUSY events is 
peaked closer to 1 than in the UED events, as one would expect given the dependence of 
the matrix elements. Following the same procedure as before we constructed the probability 
distributions f{qH\H), and determined the likelihood ratio r for each experiment. This r 
distribution is given in Fig. |8] for both SUSY and UED events. 

Comparing Fig. |8] and Fig. |4] one can see that by combining both cos and (3j, it is 
possible to enhance the separation of the SUSY and UED events, compared to reconstructing 
cos 6*^ alone. Defining the acceptance region for the SUSY hypothesis as: 

ln(r) > 16. (54) 

we obtain a statistical test with significance a = 1.2% and power 1 — /3 = 99.1%. This test 
improves over using cos6'^ alone by a factor of ^ 2, and is over 10 times less likely to result 
in error than the cos 0u method. 

Attempts to further improve the hypothesis test by including the kinematic variables 0, 
9f^ and 6^ were unsuccessful. Although our method could reconstruct some features of these 
distributions they were not statistically significant once the background contribution was taken 
into account. 



3.3.1 The effect of mass uncertainties 

Errors in the measurement of the masses and will affect the efficacy of our proposed 
statistical tests. To study the effect of such a systematic error we varied the true mass 
m^, used to generate the events, whilst keeping our statistical test fixed. (That is, the 
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Variable 


Acceptance Region 


m_4 = m_4 


rriA = 0.9m^ 


m_4 = 1.1m A 


a 


1-/3 


a 


1-/3 


a 


1-/3 


qn {cos 0a) 
qnicosOA, f^A) 
cos 6ii 


ln(r) > 2 
ln(r) > 16 
\n{rii) > 


3.3% 
1.2% 
14.8% 


98.3% 
99.1% 
84.9% 


3.3% 
2.8% 
12.3% 


96.1% 
99.7% 
84.6% 


4.7% 
0.5% 
15.3% 


98.5% 
95.2% 
82.7% 



Table 2: The significance a and power 1 — /? of several statistical tests used to discriminate between SUSY 
and UED models, for various true masses niA- 



masses assumed for the reconstruction were fixed at = 50 GeV and = 350 GeV, 
and we used the same probability densities f{qH\H) and acceptance region defined in the 
previous section.) Assuming an 0(10%) uncertainty in the mass m^, we performed another 
2000 pseudo-experiments for the cases = IAttia and riiA — 0.9mA- For each case, the 
significance and power of each statistical test is given in Table [2] 

We observe that the tests based upon the qnicosOA) and cos On distributions are quite in- 
sensitive to the 10% mass uncertainty, with variations in a and /3 at the level of 10 — 50%. Thus 
these tests, at least, are fairly robust against uncertainties in the mass measurement. How- 
ever, the test based upon qnicosOA, 13a) exhibits a much larger variation with /3 = 0.9^0 6% 
and a = 1.2~^q'^%. This large variation appears to be caused by the sensivity of the QhWa) 
distribution to the true and assumed masses. Hence, a 10% mass uncertainty reduces the 
significance and power of this test, causing it to be less effective than the test based upon 
9/f(cos 9 a)- 

4 Conclusions 

In this paper we have proposed a new method for reconstructing the probability distributions 
of kinematic variables that are sensitive to the spin of new particle states. This method relies 
on performing a probabilistic reconstruction of each event, combining information from the 
measured particle masses and the matrix element for a hypothesised process. As a result, one 
is able to study kinematic observables that depend on unknown momenta, such as that carried 
by two potential DM candidates. By combining the information from a sample of events we 
have shown how one can attempt to reconstruct this probability distribution and then test 
the assumed hypothesis. 

We have performed a preliminary test of the method on a candidate process: the discrim- 
ination of SUSY and UED models in slepton-pair production. This test was performed at 
parton level, without hadronisation or detector simulation, and took into account the effects 
of event selection and the dominant SM backgrounds. We have shown that with 100 fb^^ 
of luminosity one can discriminate between the SUSY and UED models, with a probability 
of error of 0(3%). This result has been shown to be robust against a 10% uncertainty in 
the mass measurement. Our complete results are summarised in Tableland show that, in 
principle, one can improve the test by combining several independent variables. However, the 
improvement gained in practice will depend upon the uncertainty in the masses. 

We compared this result with that expected for an alternative method, which uses the 
variable cos On. We find, for the process considered, that this can perform a discrimination 
with a probability of error of 0(15%). Thus, potentially, our method offers an improved way 
of discriminating between models and measuring the spin. 

Despite this encouraging result, many further questions remain and deserve detailed study. 
The true efficacy of our method can only be understood through a real world application, 
together with more detailed simulation. In particular, it will be crucial to study the effect of 
systematic uncertainties in the background estimation and mass measurement, combined with 
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the effects of detector resolution and acceptance. Furthermore, it would be of interest to study 
how this method performs in a wider range of processes. Possible processes of interest are 
those that involve long decay chains, or those involving diagrams with both s- and t-channel 
particle exchange. 
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